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Abstract. In the present paper the three state Potts model with competing binary interactions (with 
couplings J and J v ) on the second order Bethe lattice is considered. The recurrent equations for 
the partition functions are derived. When J v = 0, by means of a construction of a special class of 
limiting Gibbs measures, it is shown how these equations are related with the surface energy of the 
Hamiltonian. This relation reduces the problem of describing the limit Gibbs measures to find of 
solutions of a nonlinear functional equation. Moreover, the set of ground states of the one-level model 
is completely described. Using this fact, one finds Gibbs measures (pure phases) associated with the 
translation-invariant ground states. The critical temperature is exactly found and the phase diagram is 
presented. The free energies corresponding to translations-invariant Gibbs measures are found. Certain 
^ — , physical quantities are calculated as well. 
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1. Introduction 



The Potts models describe a special and easily denned class of statistical mechanics models. Never- 
theless, they are richly structured enough to illustrate almost every conceivable nuance of the subject. 
In particular, they are at the center of the most recent explosion of interest generated by the confluence 
of conformal field theory,percolation theory, knot theory, quantum groups and integrable systems. The 
Potts model |Poj was introduced as a generalization of the Ising model to more than two components. 
At present the Potts model encompasses a number of problems in statistical physics (see, e.g. |W|). 
Some exact results about certain properties of the model were known, but more of them are based on 
approximation methods. Note that there does not exist analytical solutions on standard lattices. But 
investigations of phase transitions of spin models on hierarchical lattices showed that they make the 
exact calculation of various physical quantities _DGM| , |P11 IP2j , jTj . Such studies on the hierarchical 
lattices begun with development of the Migdal-Kadanoff renormalization group method where the lat- 
tices emerged as approximants of the ordinary crystal ones. On the other hand, the study of exactly 
solved models deserves some general interest in statistical mechanics |Baj . Moreover, nowadays the 
investigations of statistical mechanics on non-amenable graphs is a modern growing topic (|L"]). F° r 
example, Bethe lattices are most simple hierarchical lattices with non-amenable graph structure. This 
means that the ratio of the number of boundary sites to the number of interior sites of the Bethe 
lattice tends to a nonzero constant in the thermodynamic limit of a large system, i.e. the ratio W n /V n 
(see for the definitions Sec. 2) tends to (k — l)/(k + 1) as n — > co, here k is the order of the lattice. 
Nevertheless, that the Bethe lattice is not a realistic lattice, however, its amazing topology makes the 
exact calculation of various quantities possible [E]. It is believed that several among its interesting 
thermal properties could persist for regular lattices, for which the exact calculation is far intractable. 
In |PLM1| IPLM2j the phase diagrams of the g-state Potts models on the Bethe lattices were studied 
and the pure phases of the the ferromagnetic Potts model were found. In jH] using those results, un- 
countable number of the pure phase of the 3-state Potts model were constructed. These investigations 
were based on a measure-theoretic approach developed in |Ge| , |Prj , |S] , |P1| IP2j . The Bethe lattices 
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were fruitfully used to have a deeper insight into the behavior of the Potts models. The structure of 
the Gibbs measures of the Potts models has been investigated in [UJ IGRj . Certain algebraic properties 
of the Gibbs measures associated with the model have been considered in [M] . 

It is known that the Ising model with competing interactions was originally considered by Elliot [Ej 
in order to describe modulated structures in rare-earth systems. In jBBJ the interest to the model was 
renewed and studied by means of an iteration procedure. The Ising type models on the Bethe lattices 
with competing interactions appeared in a pioneering work Vannimenus [V], in which the physical 
motivations for the urgency of the study such models were presented. In |YOSl ITYj the infinite- 
coordination limit of the model introduced by Vannimenus was considered. It was also found a phase 
diagram which was similar to that model studied in BB . In |MTA| . [SU] other generalizations of the 
model were studied. In all of those works the phase diagrams of such models were found numerically, 
so there were not exact solutions of the phase transition problem. Note that the ordinary Ising model 
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jBTTl 


BRZ1 


BRZ2, 


BRSSZ 



investigated. In |GPW11 lGPW2j . |MRll |MR2| the Ising model with competing interactions has been 
rigourously studied, namely for this model a phase transition problem was exactly solved and a critical 
curve was found as well. For such a model it was shown that a phase transition occurs for the medium 
temperature values, which essentially differs from the well-known results for the ordinary Ising model, 
in which a phase transition occurs at low temperature. Moreover, the structure of the set of periodic 
Gibbs measures was described. While studying such models the appearance of nontrivial magnetic 
orderings were discovered. 

Since the Ising model corresponds to the two-state Potts model, therefore it is naturally to consider 
g-state Potts model with competing interactions on the Bethe lattices. Note that such kind of models 
were studied in |NSj . |Ma] . |Mol| IMo2j on standard Z d and other lattices. In the present paper we 
are going to study a phase transition problem for the three-state ferromagnetic Potts model with 
competing interactions on a Bethe lattice of order two. In this paper we will use a measure-theoretic 
approach developed in |Ge| Ej, which enables us to solve exactly such a model. 

The paper is organized as follows. In section 2 we give some preliminary definitions of the model 
with competing ternary (with couplings J and J p ) and binary interactions on a Bethe lattice. In 
section 3 we derive recurrent equations for the partition functions. To show how the derived recurrent 
equations are related with the surface energy of the Hamiltonian, we give a construction of a special 
class of limiting Gibbs measures for the model at J p = 0. Moreover, the problem of describing the limit 
Gibbs measures is reduced to a problem of solving a nonlinear functional equation. In section 4 the 
set of ground states of the model is completely described. Using this fact and the recurrent equations, 
in section 5, one finds Gibbs measures (pure phases) associated with the translation-invariant ground 
states. A curve of the critical temperature is exactly found, under one there occurs a phase transition. 
In section 6, we prove the existence of the free energy. The free energy of the translations-invariant 
Gibbs measures is also calculated. Some physical quantities are computed as well. Discussions of the 
results are given in the last section. 



2. Preliminaries 

Recall that the Bethe lattice T k of order k > 1 is an infinite tree, i.e., a graph without cycles, such 
that from each vertex of which issues exactly k + 1 edges. Let T k = (V, A), where V is the set of 
vertices of F k , A is the set of edges of T k . Two vertices x and y are called nearest neighbors if there 
exists an edge I G A connecting them, which is denoted by / =< x,y >. A collection of the pairs 
< x, x\ >,...,< Xd-i,y > is called a path from x to y. Then the distance d(x,y),x,y £ V, on the 
Bethe lattice, is the number of edges in the shortest path from x to y. 
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Figure 1 . The first levels of T\ 

For a fixed x° £ V we set 

W n = {x £ V\d(x, x°) = n}, V n = U™ =1 W m , 
Ln = {I =< x, y >£ L\x, y £ V n }. 

Denote 

S(x) = {y£ W n+ i : d{x, y) = 1}, x £ W n , 

this set is called a set of direct successors of x. 

For the sake of simplicity we put |x| = d(x, x°), x £ V. Two vertices x, y £ V are called the second 
neighbors if d(x,y) = 2. Two vertices x,y £ V are called one level next-nearest-neighbor vertices if 
there is a vertex z £ V such that x, y £ S(z), and they are denoted by > x,y <. In this case the 
vertices x,z,y are called ternary and denoted by < x,z,y >. In fact, if x and y are one level next- 
nearest-neighbor vertices, then they are the second neighbors with \x\ = \y\. Therefore, we say that 
two second neighbor vertices x and y are prolonged vertices if |x| ^ \y\ and denote them by > x,y <. 

In the sequel we will consider semi-infinite Bethe lattice of order 2, i.e. an infinite graph without 
cycles with 3 edges issuing from each vertex except for x° that has only 2 edges. 

Now we are going to introduce a semigroup structure in (see |FNWj ). Every vertex x (except 
for x°) of has coordinates (ii, . . . , i n ), here i^ G {1, 2}, 1 < k < n and for the vertex x° we put (0). 
Namely, the symbol (0) constitutes level and the sites (ix, . . . ,i n ) form level n of the lattice, i.e. for 
x £ r^_, x = (ii, . . . ,i n ) we have \x\ = n (see Fig. 

Let us define on T^. a binary operation o : x — > as follows: for any two elements 
x = (h, . . . , i n ) and y = (j 1 ,... ,j m ) put 

(2.1) xoy = (i 1; . ..,i n )o . . ,j m ) = (i x , . . .,i n ,ji,. . . ,j m ) 
and 

(2.2) x o x° = x° o x = (h, ...,i n )o (0) = (h, . . .,i n ). 

By means of the defined operation becomes a noncommutative semigroup with a unit. Using 
this semigroup structure one defines translations r g : T 2 + — > T 2 + , g £ T\ by 

(2.3) T g (x)=gox. 
It is clear that T( ) = id. 

Let 7 be a permutation of {1, 2}. Define n^l} : r 2 + - T% by 



(2.4) 



'(0) 

7T ( ( >) = (0) 

■K { $(ii,...,i n ) = (7(»l] 
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for all n > 1. For any g G (5 7^ x°) define a rotation 7Tg 7 ^ : — ► T 2 ^ by 

(2.5) vrW(^)=r g (vr ( ( J ) (x)), x G r+. 

Let G C be a sub-semigroup of and /i : — > R be a function denned on r+. We say 
that /i is G-periodic if h(r g (x)) = h(x) for all g G G and x G Any r+-periodic function is called 

translation invariant. We say that h is quasi G-periodic if for every g G G one holds /i(vrj 7) (x)) = /i(x) 
for all x G except for a finite number of elements of T 2 ^. 
Put 

(2.6) G fe = {xGr^: |x|/fcGN}, k>2 

One can check that G^ is a sub-semigroup with a unit. 

Let <& = {771,772, ...,77^}, where 771, 772, 77^ are elements of R g_1 such that 



(2-7) 77^ 



1, for i = j, 



7-1 



here xy, x,y G M 9 , stands for the ordinary scalar product on 
From the last equality we infer that 

g 

(2.8) J>fc = 0. 

k=i 

The vectors {771, 772, 775-1} are linearly independent, therefore further they will be considered as 
a basis of W~ l . 

In this paper we restrict ourselves to the case q = 3. Then every vector /1GK 2 can be represented 
as h = h\?]i + /12772, i.e. h = (hi, /12), and from Q2.7J) we find 

{hi - \h%, if i = l, 

-\h x + h 2 , if i = 2, 

-\{hi + h 2 ), if i = 3. 

Let = (V, A). We consider models where the spin takes its values in the set <3? = {771,772,773} 
and is assigned to the vertices of the lattice T 2 ^. A configuration a on V is then defined as a function 
x G V — ► cr(x) £ in a similar fashion one defines configurations a n and cx^ on and W n , 
respectively. The set of all configurations on V (resp. V n , W n ) coincides with Q = <fr v (resp. VLy n = 
$ n , Q\y n = ). One can see that 0y„ = 0y„_i X 0w/„- Using this, for given configurations 
c n _i G 0y n _i and cj( n ) G &w„ we define their concatenations by the formula 

^n-iVaW = !{a n (x),x G K-i}, {^ (n) (y), V € VF„} 



It is clear that cr n _i V cj( n ) G 0y„. 

The Hamiltonian of the Potts model with competing interactions has the form 



( 2 -10) H(a) = —J' 8 a {x)(T(y) — Jp X/ ^{x)a{y) - J[ ^ ^(a)o-(w) 

>z,S/< <*,3/> 

where J', J p , J{ gl are coupling constants, a G and 5 is the Kronecker symbol. 
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3. The RECURRENT equations for the partition functions and Gibbs measures 

There are several approaches to derive an equation describing the limiting Gibbs measures for the 
models on the Bethe lattices. One approach is based on properties of Markov random fields, and 
second one is based on recurrent equations for the partition functions. 

Recall that the total energy of a configuration a n 6 0y n under condition a n G &v\V n is defined by 

H(a n \a n ) = H(a n ) + U(a n \a n ), 

here 

H(a n ) = -J' 6 <Tn(x)<J n (y) - Jp Y 6 <rn(.x)a n (y) 

> x,y < >~x~y~< 

x >y e v n x,y ev n 



- J 'l Y S °n(x)>Tn 



(y) 



< x, y > 

x , y G V„ 



U(a n \a n ) = -J' Y S <T n (x)* n (y) ~ J P Y S < 



> x,y < > x,y < 

x e v n , x e v n , 

v ev\v n y ev\v n 



( 3 - 2 ) - J l Y S <Tn(x)a r 



(y) 



< x, y > 

x e v„, 
v e v \ v n 



The partition function in volume V n under the boundary condition a n is defined by 

Z^= exp(-[3H(a\a n )), 

ae0 Vn 

where j3 = 1/T is the inverse temperature. Then the conditional Gibbs measure /j„ in volume V n 
under the boundary condition a n is defined by 

fin{<r\<r n ) = ^ , cr e Vn . 

Consider 0y 1 - the set of all configurations on V\ = {(0), (1), (2)}, and enumerate all elements of it 
as shown below: 

^ 9{ *~ 1)+1 = {m,m,Vi}, cr 9 ^- 1 ^ 2 = {^,771,772}, cr 9 ( i - 1 )+ 3 = {ru, m,Vs}, 

(7 9( l -l)+4 = {77^,772,171}, CT 9 (*- 1 )+ 5 = {77j, T/2, 772}, CT 9 ( ?_1 ) +6 = {f]i, 772, 773}, 
<T 9(i-l)+7 = {77^ 773, 7/i}, ,j9(i — 1)+8 = r?3) r?2 } ) cr9i = {77^, 773, 7/ 3 }, 

where i G {1, 2, 3}. 

We decompose the partition function Z n into 27 sums 

i=i 

where 

= Y exp(-PH n (a n \a n )), i € {1, 2, ■ ■ ■ , 27}. 

We set 

= exp(/3J'); # P = exp(/3J p ); ^ = exp(/3J(); 
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and 
(3.3) 

that is 



exjp(-f3H n (a n \cr n )), ie {1,2,3}, 



z 



(») 

9(i-l)+fc> 



ie {1,2,3}. 



Taking into account the denotation (|A.1I) through a direct calculation one gets the following system 
of recurrent equations 



(3.4) 





H) 


= ^?(4 n) ) 2 > 


7 (n+l) 
z 10 


= ^(4 n) ) 2 , 


7 (n- 
Z 19 


-1) 


= *(4 B) ) 2 . 


7 {n- 
Z 2 


H) 


= 1 4 n) J B 1 n) , 


7 (n+l) 
z ll 




7 (n- 
Z 20 


hi) 


= 4 n) B 3 >) , 


Z 3 


W) 




7 (n+l) 
z 12 


— jiWrW 

— ^2 °2 ' 


7 (n- 
Z 21 


HI) 


= fliA( n) 4 n) 


7 {n+l) 
Z 4 




7 (n+l) 
z 13 




7 (n- 
Z 22 


Hi) 


_ 7 ("+i) 

— ^20 ' 


7 (n- 
Z 5 


M) 


= ^(V } ) 2 , 


7 (n+l) 
Z 14 


= ee\{B^f, 


7 (n- 
Z 23 


HI) 


= o(B^r, 


Z 6 


H) 




7 (n+l) 
Z 15 




7 (n- 
Z 23 


HI) 






W) 


_ 7 (n+l) 
— z 3 ' 


7 (n+l) 
z 16 


— Z 12 > 


7 (n- 
Z 25 


HI) 


_ 7 0+i) 

— Z 21 > 


Z 8 


H) 


_ 7 (n+l) 


7 (n+l) 
Z 17 


— Z 15 ' 


7 (n- 
Z 26 


HI) 


_ y(n+l) 
— ^24 > 


Z 9 


H) 


= ^(Ci n) ) 2 , 


7 (n+l) 
Z 18 


= ^(^ n) ) 2 , 


7 (n- 
Z 27 


HI) 





Introducing new variables 



(3.5) 



(n) 
1 


= Z[ n) ; 


(n) 

x 2 


_ 7 (n) 
- Z 2 


- ^4 > 


(n) 
x 3 


- Z 3 


(n) 
4 


- Z 5 ' 


(n) 
x 5 


- Z 6 


_ 7 ( n ). 


(ra) 

x 6 


- Z 9 


(») 
7 


— Z 10 i 


(n) 

x 8 


- Z ll 


- Z 13 ' 


(n) 
X g 


- z 12 


(n) 
10 


_ y( n ). 
— Z 14 ! 


(n) 
x ll 


- Z 15 


_ y( n ). 
- Z 17 ' 


(n) 
x 12 


_ y( n ) 
— z 18 


(») 
13 


_ 7 ( n ). 
— Z 19 ! 


(n) 

x M 


_ y( n ) 
— Z 20 


_ y( n ). 
— Z 22 ) 


(n) 
x 15 


_ y( n ) 
— Z 21 


(n) 
16 


_ y( n ). 

— Z 23 ' 


(n) 
X 17 


_ 7 ( n ) 
— Z 24 


_ 7 ( n ). 
— Z 26 ' 


(n) 
x 18 


_ y( n ) 
— Z 27 



(«). 
16 ' 



y 25 » 



the equations (|3.4|) are represented by 



(3.6) 



r r (n+l) 

x l 


= eel(Af ] ) 2 


x (n+l) 


= OrA^B^ 


T (n+1) 
3 


= e 1 A[ n) c ( { l \ 


Jn+1) 


= 0{Bt ] )\ 


x 5 


= B[ n) c[ n \ 


x 6 


= e(c{ n) r, 




= 0(4 n) ) 2 , 


T (»+l) 
8 


= 9xA^B^ 


< T ("+ 1 ) 
9 


_ jW^W 

— ^2 L/ 2 ) 


r ("+!) 
x 10 


= eeliB^f 


x ll 


= e 1 B^ct\ 


T (»+l) 
x 12 


= e(ct ] ? 


T ("+l) 
x 13 


= ^(4 n) ) 2 , 


r ("+l) 
x 14 


— /l( n ) rW 

— ^3 -°3 ; 


x 15 


= 9 1 Al l) ct ) , 


r (»+ 1 ) 
x 16 


= 0{B^)\ 


T (»+l) 
v x 17 


= e 1 B^c^ n) , 


T ("+l) 
x 18 


= 00?(cf } ) 2 



A- 



The asymptotic behavior of the recurrence system ()3.6|) is defined by the first date {x2"' 
1, 2, . . . , 18}, which is in turn determined by a boundary condition a. 

Let us separately consider free boundary condition, that is U(a\a) is zero, and three boundary 
conditions a n = r]i, where i = 1,2,3. Here by a n = rj we have meant a configuration defined by 
a n = {a(x) : a(x) = rj, Vx G V \ V n }. 
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For the free boundary we have 

(l) _ a . 

(1) _ 

6 — 

(1) _ 

9 — 

(1) _ 

12 ~~ 

(1) _ 

15 — 

(1) _ 

18 — 



Jl) 

x l 


= 88 2 ; 




= 9r, 


X 


X4 


= 8; 


r h 

x 5 


= i; 


X 


T (l) 


= 8; 


T (l) 


= 9r, 


X 


T (l) 
x 10 


= 891; 


T (l) 

x ll 


= 9r, 


X 


r (l) 
x 13 


= 0; 


T (l) 

X 14 


= i; 


X 


T (l) 
x 16 


= 9; 


x 17 


= 9r, 


X 



and from the direct calculations (see (|A.2|l ) we infer that 



(n) 
1 

(n) 



B. 



(n) 
2 

(n) 



c. 



AT' = B, 



(n) 
3 ' 
in) 



(n) 



a 



(n) 



so that 



y(n) _ y(n) _ y(n) 

Zj -1 Zjry Zj O 



Jl — tj 2 — ^3 

Hence the corresponding Gibbs measure (jlq is the unordered phase, i.e. /j,(a(x) = 7/j) 
xeT%,i = 1,2,3. 

Now consider boundary condition a = rjx. Then we have 



1/3 for any 



(1) 

1 


qq6q4. 

— vV x V p , 


2 


- 9f9p] 


T (i) 

3 


(1) 




r w 

x 5 




r (i) 

x 6 


4 


= 99p, 


= 9p, 


(1) 
7 


= 99f; 


T (l) 
x 8 


= 9f; 


r (l) 
9 


(1) 




r (l) 
x ll 




r (l) 
x 12 


10 


= 99f; 


= 9t; 


(1) 




T (l) 
x 14 




r (l) 
x 15 


13 


= 06$; 


= 91; 


(1) 




T (l) 
x 17 




T (l) 
x 18 


16 


= 9; 


= 9t; 



9\9p; 

89 

9b 



8h 

eel 



By simple calculations (see (|A.2|0 we obtain 



B 
B. 



(ft) 
1 

in) 



(ft) 



I (ft) 

l 2 

(n) 



C, 



3 > 
(ft) 



and 



(ft+1) 



^(A[ n) ) 2 + 40^"^ + 2(9 + 1)(B\ 



1 (ft) d(") 



W\2 



(ft+1) 



z. 



(ft+1) 



26 x Af ) Bf ) + 2Af ) C% l) 



e(A^f 

+99 2 {B^ ) ) 2 + 2B X B^C^ + #(C^ n) ) 2 - 



By the same argument for the boundary condition a = 772 we have 



(ft) 



Kft) 



and for the boundary condition (7 = 7/3 



If 6* p = 1, i.e. J p 



y i — ^2 • 

0, then from the system of equations (|3.4|) we derive 



Z 



(ft+i) 
1 



(3.7) Z (n+1) 



00 2 (Z) 



{n) ^ 2 + 20 1 Z 1 (n) ^ n) + 26 x Z ( fzf ) + fl(Z< n) ) 2 + 2^^ n) ^ n) + 9(Z^ 2 



2 

7-(ft+l) 
Z 3 



;(n) y(n) 



r(n) y(n) 



9{Z { ? ) ) 2 + 29 1 Z ( / l) zi n) + 2z{ n) Z 3 n) + 99 2 (zf ) ) 2 + 20^" } z] n) + 0(Z (n) ^ 2 
9(z[ n) ) 2 + 2z[ n) Z^ n) + 20^ Z^ + 9{Z '" 



j-(ft) J7( n ) 
Z 3 



32^ ^( n h2 

J 2 



3 

(ft)\2 
3 ) 
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Letting 



then from (|3.7j) one gets 



Z 



(n) 



u,. 



■An) 



and v r , 



Z 3 



(3. 



_ 60lul+28 1 u n v n +0vl+28 1 u n +2v n +8 
- 1 ~ 6ug+2u n v n +6vl+26 1 u n +26 1 Vn+00j ' 



Vfl+1 



8ul+20 1 u n v n +00lvl+2u n +20 1 v n +0 
0ul+2u n v„+0vl+20 1 u n +28 1 v n +00'( ' 



From the above made statements we conclude that 

(i) u n = v n = 1, Vn G N for the free boundary condition; 

(ii) v n = 1, Vn G N for the boundary condition a = 771; 

(iii) it n = 1, Vn G N for the boundary condition = r}2\ 

(iv) n n = u n , Vn G N for the boundary condition a = r/3 . 

Consequently, when J p = we can receive an exact solution. In the next section we will find an 
exact critical curve and the free energy for this case. 

Now let us assume that J p 7^ and a = 771. Then the system (|3.6(l reduces to a system consist- 
ing of five independent variables (see Appendix A), but a new recurrence system still remains rather 
complicated . Therefore, it is natural to begin our investigation with the case J p = 0. In the case 
J p 7^ a full analysis of such a system will be a theme of our next investigations |GMMPj . where the 
modulated phases and Lifshitz points will be discussed. 

Now we are going to show how the equations (|3.8[) are related with the surface energy ()4.3|) of the 
given Hamiltonian. To do it, we give a construction of a special class of limiting Gibbs measures for 
the model when J p = 0. 

Let us note that the equality (|2.7j) implies that 

2 ( 1 

Sa(x)a(y) = 3 ( <r(x)<r(y) + - 

for all x, y G V. Therefore, the Hamiltonian H(a) is rewritten by 



(3.9) 



2 2 
where J = - J , Ji = - Ji . 

3 3 1 



H(a) 



-J 



E 

>x,y< 



a(x)a{y) - Ji 



E 

<x,y> 



a(x)a(y), 



Let h : x — > h x = (h\^ x , fi2 :X ) G M 2 be a real vector-valued function of x G V. Given n 
consider the probability measure u( n ) on <& Vn defined by 



1,2,... 



(3.10) 
where 



^(a n ) = (Z^)- 1 exp{-(3H(a n )+ £ h x a n {x)}, 

xew n 



E 

>x,y<:x,y£V„ 



H{a n ) = -J 

and as before f3 = ^ and <7 n G 0y n and i?( n ) is the corresponding partition function: 



a n (x)a n {y) - Ji ^ a n (x)a n (y), 

<x,y>:x,y£V n 



(3.11) 



Z (n)= Z W(^ jfc ) = £ exp{-/3i/(<T n ) + J] ^a n (x)}. 
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Let V% C V2 C ... |J V n = V and fj,^ \ \ ... be a sequence of probability measures on $ Vl , $ 2 , ... 

n=l 

given by (|3.1U|) . If these measures satisfy the consistency condition 
(3.12) ^^K-iV^^^^K- 



-V: 

r(») 



where = {o~(x), x £ W n }, then according to the Kolmogorov theorem, (see, e.g. Ref. jHEl) there is 
a unique limiting Gibbs measure /i on (0,F), where T is a cr-algebra generated by cylindrical subset 
of 0, such that for every n = 1, 2, ... and a n G n the following equality holds 



One can see that the consistency condition (|3.12fl is satisfied if and only if the function h satisfies 
the following equation 



(3.13) 



h r X}1 = lo g F(h' y ,h' z ) 

h' xa = \o g F((h' y )\{h> z n 

here and below for given vector h = (hi, /12) by h! and h f we have denoted the vectors |/i and (J12, hi) 
respectively, and F : R 9_1 x R 9_1 — > R is a function defined by 

]20 e hi+n + Oi(e hl+r2 + e h2+ri ) + 9e h2+r2 + 6q(e hl + e ri ) + e h2 + e ra + 9 



(3.14) F(/i,r) 



where h = (hi, h<i),r = (ri,T2) and < y,x, z > are ternary neighbors (see Appendix B for the proof). 

Consequently, the problem of describing the Gibbs measures is reduced to the description of solutions 
of the functional equation ()3.13|) . On the other hand, we see that from the derived equation (|3.13|) we 
can obtain (|3.8|) . when the function h is translation invariant. 



4. Ground states of the model 

In this section we are going to describe ground states of the model. Recall that a relative Hamiltonian 
H (a, ip) is defined by the difference between the energies of configurations a, ip 

H(a,(f) = - J' ^ (S a (x)a{y) - $<p(x)<p(y)) ~ A ^ ^(r{x)(r{y) ~ $ <p{x) V {y)) ■> 
>x,y< <x,y> 

where J = (J 1 , J[) G R 2 is an arbitrary fixed parameter. 

In the sequel as usual we denote the cardinality number of a set A by \A\. A set c consisting of 
three vertices {xi,{x2,x^}} is called a cell if these vertices are < X2,xi,x^ > ternary. In this case, 
the vertex xi is called the origin of a cell c. By C the set of all cells is denoted. We say that two c 
and d cells are nearest neighbor if |c n d\ = 1, and denote them by < c, d >. From this definition we 
see that if c and d cells are not nearest neighbor then either they coincide or disjoint. Let a E and 
c G C, then the restriction of a configuration a to c is denoted by a c , and we will use to write elements 
of a c as follows 

a c = {cr(xi),{cr(a;2),cr(x3)}}. 

The set of all configurations on c is denoted by C . 

The energy of a cell c at a configuration a is defined by 



(4.2) 



U(a c ) - -J' K(x)a(y) ~ J'l ^ ^(x)a(y)- 

>x,y<:x,y£c <x,y>:x,y£c 
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From Q4.2J) one can deduce that for any c G C and a G we have 

C7(a c ) G {C/i(J),C/ 2 (J),C/3(J),C4(J)}, 

where 

(4.3) Ui(J) = —2J[ — J', U 2 {J)=-J' X , U 3 (J) = -J', U 4 (J) = 0, J=(J',J[). 

Denote 

Bi = {a c e0 c :U(a c ) = Ui}, i = 1,2,3,4, 
then using a combinatorial calculation one can show the following 



(4.4) 


Bl = \ 




(4.5) 


B 2 = \ 


^{vi,{r]i,Vj}},{vi,{vj,Vi}}, i^j, y'e {1,2,3} 


(4.6) 


B 3 = \ 


, Ofr, {%»?;}}> e {l,2,3}|, 


(4.7) 


Ba = \ 


|{ f 7i>{ 7 7j> r ?fc}}> G {1,2,3}, i-j-k = 6^. 


From P~T1) 


we infer that 




(4.8) 




H{^a)=Y J (U(vc)-U(o c )). 



Recall (see [R]) that a configuration G is called a ground state for the relative Hamiltonian of 
Fif 

(4.9) U(ip c ) = min{*7i(J), f/ 2 (J), f/ 3 (J), t/ 4 (J)}, for any c G C. 

A couple of configurations a,(p G coincide almost everywhere, if they are different except for a 
finite number of positions and which are denoted by a = ip [a.s]. 

Proposition 4.1. A configuration ip is a ground state for H if and only if the following inequality 
holds 

(4.10) H(ip,a)<0 
for every a £ with a = ip [a.s]. 

Proof. The almost every coincidence of a and ip implies that there exists a finite subset L C V such 

oo 

that o~(x) ^ ip(x) for all x G L. Denote Vl = f] \Yk '■ L G Vk}- Then taking into account that ip is a 

k=l 

ground state we have U{<p c ) < U(a c ) for every c G C. So, using the last inequality and (|4.8|) one gets 

H(ip,a)= (U(<p e )-U(v c ))<0. 

c&C,c&V L 

Now assume that (|4.1Uj) holds. Take any cell c G C. Consider the following configuration: 

n ( T \-{ (T ( X )' if x G c > 
^c )V ^J - | ^ x)) if ^ ^ Cj 

where a G C . It is clear that a C)tp = ip [a.s.], so from (|4,8[) and (j4,10|) we infer that H(ip,a Cttp ) = 
U(ip c ) — U(o~) < 0, i.e. U(ip c ) < U(a). From the arbitrariness of a one finds that p is a ground 
state. □ 
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Denote 

A k = {jGM 2 : U k (J)=mm{U l (J),U 2 (J),U 3 (J),U A (J)}}, k = 1,2,3,4. 

From equalities (|4.3|) we can easily get the following 

A 1 = {J = (J', J[) G R 2 : J[ > 0, J[ + J'> 0} 

A 2 = {J = (J', J[) G R 2 : J[ > 0, J[ + J' < 0} 

,43 = { J = (J', J{) g M 2 : J[ < 0, J' > 0} 

A 4 = {J = (J', J{) G M 2 : J[ < 0, J' < 0} 



Denote 



A k \(\jA k n Aj\ 
S'=i ' 



B k = A k \[ \A k nAj), k = 1,2,3,4. 



Now we are are going to construct the ground states for the model. Before doing it let us introduce 
some notions. Take two nearest neighbor cells c, d G C with common vertex x G cf)c'. We say that two 
configurations a c G C and a c > G C / are consistent if cr c (x) = <r c /(2;). It is easy to see that the set V can 
be represented as a union of all nearest neighbor cells, therefore to define a configuration a on whole 
V, it is enough to determine one on nearest neighbor cells such that its values should be consistent on 
such cells. Namely, each configuration a G is represented as a family of consistent configurations 
on C , i.e. a = {<7 c } ce c- Therefore, from the definition of the ground state and (|4.4j) - ()4,7[) we are able 
to formulate the following 

Proposition 4.2. Let J G B k then a configuration if = {f c }ceC is a ground state if and only if 
if c G B k for all c G C. 

Let us denote 

(4.11) <7^ m ) = {o~(x) : o~(x) = r] m , \/x G V}, m = 1,2,3. 

Theorem 4.3. Let J G Bi, then for any fixed a c G Bi (here c is fixed), there exists a ground state 
ip G with ip c = a c . 

Proof. Let a c G Bi. Without loss of generality we may assume that the center x\ of c is the origin of 
the lattice T 2 ^. Further we will suppose that cr(xi) = n± (other cases are similarly proceeded). Put 



N®(c. 



c I 



3 = 1,2,3, 



j*; G {1,2,3} : a c (x k ) = rjA 
(VfVc),AfVc),AfVc)), CGC. 



It is clear that ivfVc) > and £ N^(a c ) = 3. 

k=i 

According to Proposition 14.21 to find a ground state (p G it is enough to construct a consistent 
family of ground states {(p c } c ^c- 

Consider several cases with respect to i (i G {1, 2, 3, 4}). 

Case i = 1. In this case, according to (|4.4|) we have cr c (x) = r/i for every x G c. This means that 
ni{o~c) = (3,0,0). Then the configuration o-W is the required one and it is a ground state. From ()2.3j) 
we see that a® is translation-invariant. 

Case i = 2. In this case from (|4,5j) we find that n,2(cr c ) is either (2, 0, 1) or (2, 1, 0). Let us assume 
that fi2(o~c) = (2,0, 1). Now we want to construct a ground state on nearest neighbor cells, therefore 
take c', c" G C such that < c, d >, < c, c" > and c' ^ c" . It is clear that d n c" = 0. Let X2 and X3 
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\ / \ ' \ / \ / 

\ / \ ' \ / \ / 

\ 1 \ ' \ / \ / 

\ / \ ' \ 1 \ / 

\ / \ * \ / \ / 

\ / \ ' \ / \ / 

\ 1 \ 1 \ 1 \ 1 




Figure 2. f^ 1 ' 3 ) — ground state. The coupling constants belong to B2 

be the centers of c' and c", respectively. So due to our assumption we find that either (t{x2) = rji, 
a(x^) = 773 or cr(x2) = 773, cr{x2) = r)\. Let us consider a{x2) = rji, a{x^) = 773. Then we have 
Cc = {771, {?7i, ^73}}- We are going to determine configurations (p c > S C ', £ 0c" consistent with a c 
and N\ (a) ■ (a) = 2, a = (f' c ,(p". To do it, by means of (|4.5I) . we choose configurations <p c and 
</V on c',c", respectively, as follows 

(4.12) 9V = {771, {771, 773}}, W = {773, {771,773}}- 

Hence continuing this procedure one can construct a configuration ip on V, and denote it by y^ 1 ' 3 ). 
From the construction we infer that <^( 1,3 ) satisfies the required conditions (see Fig. EJ). The constructed 
configuration is quasi T^-periodic. Indeed, from (|2.4|) and (|4.12|) one can check that for every x £ Tj_ 

with \x\ ^ 1 we have p (1,3) (?r£j)(aO) = v? (1 ' 3) (z), here 7({1,2}) = {2, 1}. So from for every 5 £ Tl 
one finds that ^(tt^^)) = v? (1 ' 3) (^) 

for all |x| 7^ 1. Similarly, we can construct the following quasi 

periodic ground states: 

^ 3 ' 1 ),^ 1 ' 2 ),^ 2 - 1 ),^ 2 ' 3 ),^ 3 ' 2 ). 

Case i = 3. In this setting we have that n 3 (cr c ) is either (1,0,2) or (1,2,0) (see ()4,6|0 . Let us 
assume that n2(<T c ) = (1,2,0). Let c',c" £ C be as above. From (j4.6|) and our assumption one 
finds a(x 2 ) = (y{xz) = r/ 2 . Then again taking into account (|4,6|) for c',c" we can define consistent 
configurations by 

(4.13) tp c > = {v2,{Vi,Vi}}, fc" = {7?i,ryi}}. 

Again continuing this procedure we obtain a configuration on V, which we denote by ^jP-' 2 l. From 
the construction we infer that <£>[ 1,2 1 is a ground state and satisfies the needed conditions (see FigEJ). 
From (|4.13jl and (|2.cS|) we immediately conclude that it is GVperiodic. Similarly, we can construct the 
following GVperiodic ground states: 

Note that on c', c" we also may determine another consistent configurations by 

(4-14) tp c > = {7?2,{773,7?3}}, fc" = {772, {7?3,7?3}}. 
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TV '^3 V ^3 V A V ,^3 

\ i \ i \ i \ i 

\ i \ 1 \ / \ / 

\ i \ ' \ / \ / 

\ / \ ' \ i \ / 

\ / \ * \ / \ / 

\ / \ ' \ / \ / 

\ i \ ' \ i \ i 

w \ ' \ / w 




Figure 3. ip^- 1 ' 3 ^ — ground state. The coupling constants belong to B% 

Now take b',b" G C such that < c',&' >, < c',b" > and b' / 6". On 6', b" we define consistent 
configurations with ip c / by 

(4.15) W = W = {V3,{Vi,Vl}}- 

Analogously, one defines 92 on the neighboring cells of c" . Consequently, continuing this procedure 
we construct a configuration ^I 1,2 ' 3 ! on V. From ([2.6j) . ([2.3|) . (|4.14|) and 1)4.15(1 we see that 

^[1,2,3] is 

a G^-periodic ground state. Similarly, reasoning one can be built the following G^-periodic ground 
states: 

^[2,1,3] ^[2,3,1] ^[1,3,2] ^[3,1,2] ^[3,2,1] . 

These constructions lead us to make a conclusion that for any number of collection . . . , i^} with 
G {1,2,3} we may construct a ground state (pfa'-^ which is G^-invariant. Hence, 
there are countable number periodic ground states. 

Case i = 4:. In this case using the same argument as in the previous cases we can construct a 
required ground state, but it would be non-periodic (see (|4.7|0 . □ 

Remark 1. From the proof of Theorem 14.31 one can see that for a given a c G Bi with i > 2, there 
exist continuum number of ground states (f G such that ip c i G Bi for any c' G C and ip c = a c . Since, 
in those cases at each step we had two possibilities there have been at least two possibilities to choice 
of ip c i and ip c n , this means that a configuration on V can be constructed by the continuum number of 
ways. 

Corollary 4.4. Let J G Bi(i ^ 4), then for any fixed a c G Bi (here c is fixed), there exists a periodic 
(quasi) ground state ip G such that ip c = a c . 

By GS(H) and GS P (H) we denote the set of all ground states and periodic ground states of the 
model (j2.1U|) . respectively. Here by periodic configuration we mean G-periodic or quasi G-periodic 
ones. 

Corollary 4.5. For the Potts model (|2,10|) the following assertions hold. 

(i) Let J G B\, then 

\GS(H)\ = \GS P (H)\=3; 

(ii) Let J G Bi then 

\GS(H)\=c, \GS P (H)\=6; 
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J' 




Figure 4. Phase diagram of ground states 

(iii) Let J E -B3 then 

\GS(H)\ = c, \GS P (H)\ = H ; 

(iv) Let J E -B4 i/ien 

|G5(fOI=c 

The proof immediately follows from Theorem 14.31 and Remark 1. 

Remark 2. From Corollary 14.51 (see Fig^J we see that when J E B\ then the model becomes 
ferromagnetic and for it there are only three translation-invariant ground states. When J £ B3 then 
the model stands antiferromagnetic and hence it has countable number of periodic ground states. The 
case J £ B2 defines dipole ground states. When J E B4 then the ground states determine certain 
solution of the tricolor problem on the Bethe lattice. All these results agree with the experimental 
ones (see |NSj ). 

5. Phase transition 

In this section we are going to describe the existence of a phase transition for the ferromagnetic 
Potts model with competing interactions. We will find a critical curve under one there exists a phase 
transition. We also construct the Gibbs measures corresponding to the ground states crW (i = 1,2,3) 
in the scheme of section 3. Recall that here by a phase transition we mean the existence of at least 
two limiting Gibbs measures (for more definitions see |Gej . jPrj ■ |S] ) ■ 

It should be noted that any transformation r g , g E r + (see (|2.3|) ) induces a shift f g : — > given 
by the formula 

(f g a)(x) = a(r g x), iGf^, a E 0. 

A Gibbs measure [i on is called translation - invariant if for every g E the equality holds 
H(f-\A)) = n{A) for all A E T, g E T 2 + . 
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According to section 3 to show the existence of the phase transition it is enough to find two different 
solutions of the equation (|3.13|) . but the analysis of solutions (|3,13|) is rather tricky. Therefore, it is 
natural to begin with translation - invariant ones, i.e. h x = h is constant for all x E V. Such kind 
of solutions will describe translation-invariant Gibbs measures. In this case the equation ()3.13|) is 
reduced to the following one 



(5.1) 



9u 2 +2uv+ev 2 +2e 1 u+2e 1 v+e'{e 

eu 2 +2e 1 uv+e 2 6v 2 +26 1 v+2u+6 
eu 2 +2uv+6v 2 +2e 1 u+2e 1 v+8' 2 6 ' 



where u = e hl ,v = e h2 for a vector h = (hi, /12). 

Thus for Op = 1 using properties of Markov random fields we get the same system of equations 

Remark 3. From (|5.1j) one can observe that the equation is invariant with respect to the lines 
u = v, u = 1 and v = 1. It is also invariant with respect to the transformation u — ► 1/u, v — > 1/v. 
Therefore, it is enough to consider the equation on the line v = 1, since other cases can be reduced to 
such a case. 

So, rewrite (|5.1|) as follows 

(5.2) u = f(u;9,9i), 

here 

,,,, ftft , fljV + 4flm + 2(0 + 1) 

(5 - 3) 05 dl) - w + 2(6i + i)u + e\e + 26i + e 

From (|5.3|) we find that (|5.2|) reduces to the following 

0u 3 + (29i - 9\9 + 2)u 2 + (0{6 + 9- 20y)u - 2(9 + 1) = 
which can be represented by 

(u - 1) ( 9u 2 + (9i + l)(9(l - 9i) + 2)u + 2(9 + 1) ) = 0. 



Thus, u = 1 is a solution of H5.2jl . but to exist a phase transition we have to find other fixed points 
of (|5.3|) . It means that we have to establish a condition when the following equation 

(5.4) 9u 2 + (9i + 1)(0(1 - 9i) + 2)u + 2(9 + 1) = 

has two positive solutions. Of course, the last one (|5.4j) has the required solutions if 

(5.5) (0i + l)(0(l-0i)+2) <0, 

(5.6) the discriminant of ()5.4|) is positive. 

The condition (|5.5j) implies that 

(5.7) 0i > 1 and 9 > 

01 — 1 

Rewrite the condition (|5.6|) as follows 

(5.8) ( (9\ - l) 2 - 8]9 2 - 4^ (0i + l) 2 (0i -I) + 2)9 + 4(0i + l) 2 > 0, 



which can be represented by 

(5.9) (e-^)(e-^)>o, 
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where 

2^(01 + l) 2 (0i-l) +2) + V(0 1 + 1)3 + 1 

From ()5.8|) we obtain that 



(5 - 10) «,;-!)* - 8 



(5.11) ft. 6- 4( " 1 + 1 



(0f_ 1)2_8' 

Now we are going to compare the condition (|5.7|) with solution of ()5.9|) . To do it, let us consider 
two cases. 

Case (a). Let (0 2 - l) 2 - 8 > 0. This is equivalent to 0i > \A + 2\/2. Hence, according to (|5.11j) 
we infer that both £1 and £2 are positive. So, the solution of (|5.9|) is 

(5.12) 0G (0,6) U (6,00). 

From 1)5.1(1 j) we can check that 

Therefore, from ()5.7|) . (|5.12j) we conclude that should satisfy the following condition 



(5.13) > 6 while 6 1 > y 1 + 2\/2. 



Case (b). Let (0 2 - l) 2 - 8 < 0, then this with (j5~T|) yields that 1 < X < Vl + 2^/2. Using (J5~9"|) 
and (|5.11|) one can find that 



(5.14) 



> £1, if 0* < 0i < Vl + 2V2 

> ^j, if 1< By < 9*, 



where 6* is a unique solution of the equation (x — l)(-\/ (x + l) 3 + 1 — 1) — 4 = 1 . 

Consequently, if one of the conditions (|5.13|) or (|5.14|) is satisfied then f(u,;9,9x) has three fixed 
points u = 1, u* and u* 2 . 

Now we are interested when both u\ and u\ solutions are attractive 2 . This occurs when 

A/( U) ; Ml )| u=1 > 1, 

du 



^One can be checked that the function 

g(x) = (x - l)(V(x + l) s + l - 1) 

is increasing if x > 1. Therefore, the equation = 4 has a unique solution 0* such that 0* > 1. 
2 Note that the Jacobian at a fixed point (u*,v*) of 15. II can be calculated as follows 



l(u*,v*) k(u*,v*) 
l(v* , u*) 



(5.15) j(u*,v*)= ( l{ f: v \ 

here 

2((0(0i -u)-(v + 0i))u + 0i (u + 1)) 



(5.16) l(u,v) 

(5.17) k(u,v) 



6u 2 + 2uv + 6v 2 + 20iu + 20i v + 0f ' 
2(1 - u)(6v + 1 + u) 

du 2 + 2uv + 6v 2 + 20i u + Wrv + e'ie ' 
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since the function f(u, ;0,0i) is increasing and bounded. Hence, a simple calculation shows that the 
last condition holds if 3 

2 



(5.18) 0i > 2 and 9 > 

9\ — 2 

If 9\ > 2 then the condition (|5.14|) is not satisfied since 9* < 2. Consequently, combining the 
conditions (j5,13|) and (|5.18|) we establish that if 

(5.19) 0i > 2 and 9 > max J 2 & 



then /(it, ;0,0i) has three fixed points, and two of them w* and 1X2 are attractive. Without loss of 
generality we may assume that u* > u 2 . Then from (|5,4j) one sees that 

u*u* - 2{ - 9 + l) 

u l u 2 — Q 

which implies that 

(5.20) u* — * 00, u 2 — > as (3 — * 00 



Let us denote 



''i.j _ (^logui,oj, h* 21 = Qlogit2,0). 



which are translation- invariant solutions of (|3.13j) . 
According to Remark 2 the vectors 



(5.21) 



h* 1>2 = (0,§loguj), h h = (0, flog u* 2 ) 

K,3 = (-fiogi4,-§io g i4), h* 23 = (-§iog«5,-|iogt4) 



are also translation-invariant solutions of (|3.13|) . The Gibbs measures corresponding these solutions 
are denoted by fiij, fJ,2,h = 1)2,3), respectively. 

From (|5.19j) we infer that (J,J\) belongs to B\. Furthermore, we assume that (|5.19|) is satisfied. 
This means in this case there are three ground states for the model. Therefore, when 8 — > 00 certain 
measures ^1^,^2,1 should tend to the ground states {cr' 1 ^, }. Let us choose those ones. Take 

then from (|XTU|) . l(23|) and lf5~2U)) we have 



e K,im + e K A ri2 + e hl 



m 



(5.22) = — J— 1 as 8 -> 00, 

v ; 1^ + 2 H 

where x £ V . 

Similarly, using the same argument we may find 

(5.23) (j,i }2 ((t(x) = 7/2) -> 1, nx^{(r{x) = 773) ->■ 1 as 3 -> 00. 

Denote these measures by = ^i,fc, = 1,2,3. The relations (|5.22|) . l|5.23|) prompt that the 
following should be true 

Hi -> (5 CT (i) as /3 -> 00, 



^Indeed, this condition also implies that the eigenvalues of the Jacobian J(l, 1) is less than one (see 15,151 - 15.171 '). 
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here 5 a is a delta-measure concentrated on a. Indeed, let us without loss of generality consider the 
measure We know that is a ground state, therefore according to Proposition 14. II one gets that 
H{ a n\v n ) > H{ a \v n ) f° r all cr G f2 and n > 0. Hence, it follows from (|3.1Uj) that 

exp{-0H(a^\ Vn ) + hl lVl \W n \} 



E exp{- i 9fr(ar n )+fc* 1 £ cr(x)} 



i+ E 



exp{-/3Jf(a n )+ftt il E *(x)} 
e«p{-^(«TCi)|v n )+A? x i7i|W„|} 



> ; ► 1 as /? — > oo. 

" 1 + 1/uJ 

The last inequality yields that the required relation. 

Consequently, the measures \ik (k = 1, 2, 3) describe pure phases of the model. 
Let us find the critical temperature. To do it, rewrite (|5.19|) as follows: 

(5 ' 24) i^AirHl 

where 

f(x) = xlog 



exp(l/x) — 2 

/ 2( (exp(l/x) + l) 2 (exp(l/x) - 1) + 2 J + V(exp(l/a;) + l) 3 + l\ 



= xlog 



(exp(2/x) -l) 2 -8 



V 



From these relations one concludes that the critical line (see Fig|S} 4 is given by 

Consequently, we can formulate the following 

Theorem 5.1. If the condition (|5.24j) is satisfied for the three state Potts model (|.3.9|) on the second 
ordered Bethe lattice, then there exists a phase transition and three pure translation-invariant phases. 

Remark 4. If we put J = to the condition then the obtained result agrees with the results 

of jPLM1llPTM2j . [H]. 

Observation. From 1)5.15(1 - 1)5. 17j) we can derive that the eigenvalues of the Jacobian at the fixed 
points (1, «*), («2)1)> (Ij^)) (( n i) _1 ; 1 )> (( n 2) _1 i (""2) X ) are reai - Therefore, in this 

case (i.e. J p = 0), there are not the modulated phases and Lifshitz points. On the other hand, the 
absolute value of the eigenvalues of the Jacobian at the fixed points (it*, 1), (1, it*) and ((-u*) -1 , (-uf) -1 ) 
are smaller than 1. The absolute value of the eigenvalues at the fixed points (it*., 1), (l,it<!>) an d 
((ita) -1 , (it^) -1 ) are bigger than 1. These show that the points (it*,l), (l,u*) and ((u*) _1 , (u*) _1 ) 
are the stable fixed points of the transformation given by 1)5.1)1 . The Gibbs measures associated with 
these points are pure phases. 

Remark 5. Recall that the a Gibbs measure [i$ corresponding to the solution h = (0, 0) is called 
unordered phase. The purity of the unordered phase was investigated in [GTT|. pR3] when J = 0. 

Note that the functions ip and £ are increasing, therefore their inverse <^ _1 and £ _1 exist. 
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Such a property relates to the reconstruction thresholds and percolation on lattices (see Mar , JM ). 
For J ^ the purity of /xo is an open problem. 

6. A FORMULA OF THE FREE ENERGY 

This section is devoted to the free energy and exact calculation of certain physical quantities. Since 
the Bethe lattice is non-amenable, so we have to prove the existence of the free energy. 

Consider the partition function Z^ n '(/3, h) (see (|3.11|) ) of the Gibbs measure //^ (which corresponds 
to solution h = {h x , x £ V} of the equation Q3.13[l 1 

zW(P,h)= Yl exp{-/?iJ(cr n ) + h x a n (x)}. 
The free energy is defined by 

(6.1) F p (h) = - lim -J—lnZ( n \(3,h). 

n-*oo Sp ■ 2 n 

The goal of this section is to prove following: 

Theorem 6.1. The free energy of the model (j3.9|) exists for all h, and is given by the formula 

1 n 

(6.2) F p (h) = - Hm zg^Yl E hga(x,hy,h g ;9,0 u p), 
where y = y(x), z = z(x) are direct successors of x; 



(6.3) a(x,h y ,h z ;0,e 1: P) = e-W 2+J ^g(ti y ,ti z ) 

here the function F(h,r) is defined as in (|3,14j) . and 

g (h, r) = 9e hl+ri + e hl+T ' 2 + e h2+Tl + 9e h2+r2 + 9^ + e ri + e h2 + e T ' 2 ) + 0\d, 
where h = (h\, r = (ri, r-i). 

Proof. We shall use the recursive equation (|B.6|) . i.e. 

Z^ =A n _ 1 Z^ n ~ 1 \ 
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where A n = \\ a(x, h y , h z ; 9, 9 X , (3) x £ V, y, z £ S(x), which is defined below. Using (jij.3j) we have 

xew n 

Thus, the recursive equation (|B.6|) has the following form 
(6.4) Z<ri(p;h)=exp( £ loga(x, h y , h z ; 9, 9 X , /3)J Z^" 1 ^, h). 

Now we prove existence of the RHS limit of Q6.2J1 . Prom the form of the function F one gets that it 
is bounded, i.e. \F(h,r)\ < M for all h,r £ M. 2 . Hence, we conclude that the solutions of the equation 
()3.13j) are bounded, i.e. \h x i\ < C for all x E V, i = 1, 2. Here C is some constant and h x = (h x ih X 2). 
Consequently the function a(x, h y , h z ; 6, Ox, /3) is bounded, and so | loga(x, hy, h z ; 0, Ox, 0)\ < Cp for 
all h y ,h z . Hence we get 



3 • 2 



1 n 

— ^2 ^2 ^ga(x,hy,h z ;9,6i,f3) 



k=i+i xew n 



(6.5) 



2 n 



k=l+l 

Therefore, from ()6.5|) we get the existence of the limit at RHS of (|6.2[) . 



□ 



Let us compute the free energy corresponding the measures /Xj, (i = 1,2,3). Assuming first that 
h x = h for all x € V. Then from (|6.2|) and (|6.3j) one gets 



^(/i) = iloga(/i,Mi,/3), 



here 
(6.6) 



a(h,e,O u P) = e-W 7 +Wg{ti i ti) 



F(h',h')F((h')\(h'y) 



1/3 



Let us consider h = h\ k , (k = 1, 2, 3). Denote = Fp[h\ k ). Then we have 



(6.7) 



= log 



-(J/2+Ji)(3(-.*\l/3(a(..*\2 



{u\yi 6 {0{ulY + 2(0! + IK + 9(9 + 20! + 



Taking into account (|5.4|) the equality (|6.8f) can be rewritten as follows: 



(6.8) pFp = -(J/2 + Jx)0 + - logu* + log(0x - 1) + log 



0(0l + l)(uJ + l) + 2 



Now let us compute the internal energy U of the model. It is known that the following formula 
holds 



(6.9) 



U 



dp ' 



Before compute it we have to calculate du\ /d/3. Taking derivation from both sides of (|5.4|) one finds 



(6.10) 



du\ 
~d/3 



3 (Ji9x( 



1) + J{9x + l))u{ + J 



29u\ + {9 l + l){9{l-9 l ) + 2) 
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From H6.8|) and H6.9|) we obtain 



U = -(J/2 + J0 + - 



Ji , 0((J + Ji)0i + J)K + i) 



+ 



1-1 

g(gi + l)(4u| + l) + 2 
_3«*(0(^i + l)(u* + l) + 2) 
Again using (|5.4j) and (|6.1U|) one gets 



i + l)« + l)+2 
d/3' 



du\ 



U 



-(J/2 + Ji) + 



ggiGMgi + 1) + - + 1) + 2J1 
-i)« + i) + 2(^-1) 



0(0i + l)(4< + l) + 2 



(0(0! + l)(00i - 2) - 2)u* - 2(0 + 1)(0! + 1) 

((j 1 e 1 (ee 1 - 1) + + i)K + j) 

20u| + (0 1 + l)(0(l-0 1 )+2) 
Using this expression we can also calculate entropy of the model. 

Since spins take values in R 2 , therefore the magnetization of the model would be M 2 -valued quantity. 
Using the result of sections 4 and 5 we can easily compute the magnetization. Let us calculate it 
with respect to the measure fi\. Note that the model is translation-invariant, therefore, we have 
Mi =< ct(o) > lll , so using (|2~9|) . (f2~8|) and (|.3.10|) one finds 



M 



i = 77i/xi(o-( ) = ??i) + ?72Mi(o-(o) = + »»Mi(^(o) = %) 



K) 2 / 3 + 2(n*)- 1 /3 



,2/3 



+ r? 2 («i 



-1/3 



+ Wi, 



-1/3 



1 



«i +2 

«J - 1 



«J + 2 



Similarly, one gets 



Mo 



M, 



< a 



(o) 



> 



1 



/'2 



itf + 2 



>72- 



1 



(0) 



2u\ + 1 
7. Discussion of results 

It is known |Baj that to exact calculations in statistical mechanics are paid attention by many of 
researchers, because those are important not only for their own interest but also for some deeper 
understanding of the critical properties of spin systems which are not obtained form approximations. 
So, those are very useful for testing the credibility and efficiency of any new method or approximation 
before it is applied to more complicated spin systems. In the present paper we have derived recurrent 
equations for the partition functions of the three state Potts model with competing interactions on a 
Bethe lattice of order two, and certain particular cases of those equations were studied. In the presence 
of the one-level competing interactions we exactly solved the ferromagnetic Potts model. The critical 
curve (|5.25|) such that there exits a phase transitions under it, was calculated (see Fig. EJ). It has 
been described the set of ground states of the model (see Fig. 0J). This shows that the ground states 
of the model are richer than the ordinary Potts model on the Bethe lattice. Using this description 
and the recurrent equations, one found the Gibbs measures associated with the translation- invariant 
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ground states. Note that such Gibbs measures determine generalized 2-step Markov chains (see |D"j K 
Moreover, we proved the existence of the free energy, and exactly calculated it for those measures. 
Besides, we have computed some other physical quantities too. The results agrees with |PLMl|fPLM 2 . 
[U] when we neglect the next nearest neighbor interactions. 

Note that for the Ising model on the Bethe lattice with in the presence of the one- level and prolonged 
competing interactions the modulated phases and Lifshitz points appear in the phase diagram (see 
fV] , |YOSj , |SCj ) . In absence of the prolonged competing interactions in the 3-state Potts model we do 
not have such kind of phases, this means one-level interactions could not affect the appearance the 
modulated phases. One can hope that the considered Potts model with J p = will describe some 
biological models. Note that the case, when the prolonged competing interaction is nontrivial (J p 7^ 0), 
will be a theme of our next investigations GMMP], where the modulated phases and Lifshitz points 
will be discussed. 
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Appendix A. Recurrent equations at J p / 



(A.l) 



4 n 

Bf 

cf 

A 2 

B f 
cf 
At 
Bf 



a 



in 



n2 7 {n) . a y(n) , Q y{n) Q y (n) y (n) y {n) „ y(n) , y(n) 



®p Z 10 



o2 y( n ) 
W p^l9 



2 -r up^ 3 -r u p ^ 4 

(«) , a y( n ) , a y{n) , q y( n ) 
' "P^ll ' "p^12 ~r "'P^13 



a y( n ) 
u P Zj 20 



a y( n ) 1 a y(n) 



6 

1 7W , y( n ) 

-t- z, 14 + z lp 

(n) 



(n) 



y\< 1 ) 1 y( n ) 1 a y( n ) _i_ y 
+ ^23 + Z 24 + °P Zj 25 + Z 26 

j2j\ + "p^2 1" ^3 ' "p^4 + "p^5 + "p-^6 7 ~"~ "p A 



+ 

An) 



P^ll 

(n) 



yV' l > 1 /) 7 V'V 1 y 
i~ t7 D Z/ 9n -f- Z/, 



P^20 



r(n) 
J 12 
r(n) , 
J 21 + 



'P' 

'P^14 
q2y( n ) 
i p Z/ 23 



+ • 



q 7 (n) , 7 (n) fl 7 (n) 



P^24 



7W 1 v( n ) 1 a y( n ) 1 7W 1 7W t a 7W , o y( n ) 



y( n ) 1 7 1 .'" ■ /i 7!"' 71"! i '7 1 "; 1 /i 7 1 ",' 1 /i 7-"'.' 1 /1 7"-; 
■^10 ' ^11 ~r t; P z '12 ~r ^13 ' z '14 ' °P^\5 T Cb^ib tt/ n ^i 7 -+- 1 



J 2 
r(n) 



r(n) 

y 12 



7W 1 7 (n) , fl 7 (n) 7 
•^19 ^20 l7 p z '9l ~r Z/r 



7 21 



4 

(n) 
13 
(n) 
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J 5 
r(n) 



J 6 

An) 



7 

(n) 
16 



/p^26 



7W 1 /1 7W 1 /1 7W , a , 
^23 + ^^24 ^^25 + ^^26 



y 17 



+ ^ 

+ Z 18 

+ ^27 

■2 y(n) 
P Z 9 
.2 7 W 

(n) 
27 



then the last one in terms of ()3.5j) is represented by 



(A.2) 



4 n 

Bf 

cf 

Af 
Bf 
Cf 



A 
B 

C 



26 p x 



(n) 



> H - *2QpXQ -\- ~h ~r" Xg 



„(«) , „(»») 



i2_(«) 
P x l 



'p-^7 

"p-^lS ' Z[7 p x 14 

xi n) + 29 p xf ] + 2xf 



+ 
+ 



(») 



("■) 1 nfl ( n ) 1 o in 

Xy + 2^ p Xs + 2x 



7 ^ * v p^8 ' ^9 
( n ) 1 9/1 ( n ) 1 9 (« 

+ 26 p xf 



(n) , m (n) 



.(«) 

» 
y 7 

.(n) 
13 



+ 2x 
+ 2x 



(n) 
2 

(n) 
8 

(n) 



+ 26 p x g 



+ 2x^ + 



r)0 ( n 
ZV P X 15 



+ 2xx;> + X 



18 ' 



+ X 



Jn) „ ( 
p-^lO ^ ^ x ll 

2^( n ) 



12 > 



p-16 
(n) 



+ 20 pX f ) + , 



-!-x 1 ( ™ ) +2^ ) + 1 



10 

+ ^16 



P 

+ 26I5QX 



'p^e 

)2™W 
p x 12 ' 



P^17 



' a p x 18 
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A{ n) 




+ (0 2 - 
' 


- l)x[ n) + 2f(9„ - 
/ 1 ^ v p 


lkl n) + 2(6„ 

/ 2 VP 




r("> 


- z (n) 

— ^2 


4- (B 2 - 
+ (Vp 


_ i -I- 9(0 — 


1 -I- 9(73 
/ x 8 + z ^"p 


IJXg 


c (n) 


- Z 3 


+ («?- 


-l) x W+2(e j ,- 


l)xff + 2(0 P 


- i)4s 




= z{ n) 


+ 2(6 P 


- i)4 n) + («? - 


l)xi n) + 2(fl p 


- i)4 n 




_ y( n ) 
- Z 2 


+ 2(6 P 


- i)4 n) + <fi - 


l)4o ) + 2(^ P 


- i)4i 




- Z 3 


+ 2(6 P 


- 1)4? + <% - 


1)4? + w P 


- l)4r 




= zj ft) 


+ 2(6 P 


- l)x^ ] + 2(9 P - 


- i)4 n) + («? 


- i)4 n 


5 w 


- Z 2 


+ 2(9 P 


- 1)4*° + 2(d p - 


- 1)4? + $ 


- 1)42 


°3 


_ v( n ) 
- Z 3 


+ 2(6 P 


- i)x { S + 2(e p - 


- 1)4? + («? 


- i)4s 



Now let us assume that J p ^ and a = rjx. Then 



B 
B. 



(n) 
1 

(n) 



(n) 



A 
B 



(n) 
2 

(n) 



A 

a 



(n) 
(n) 



and 



(n) 



Hence the recurrence system Ij3.6|) has the following form 



(A.3) 



(n+l) 



J w +1) 
J' 9 

(n+l) 



.r 



13 



(n+l) 

.(n+l) 
'17 



x (n+1) 

*(4 B) ) a , 

(n) „(n) 



x 5 
(n+l) 



,4V ; c. 



(n) „(n) 



l-°2 ^2 
(n+l) 



Jl+l) 



X 

J"+l) 
x 10 

x 12 

(n+l) 
x 14 
(n+l) 



X 



11 



X 



16 
(n+l) 
18 



x 3 

(^2 



(n) „(n) 



00f(B< 

„("+!) 
x 9 , 

Jn+1) 

x 12 ) 
„>+l) 

x 10 



Through introducing new variables 



2/1 


(n) 
= X^ 






(n) 
2/3 


(n) 
= x 5 


(n) 
3-4 


J M 
1 2/4 


(n) 


_ J") 
— x 13 


(n) 
2/6 


(n) 
— x 9 


— x 14 


(n) 
I 2/ 8 


(n) 

— x n 


_ T (») 



(n) 



(n) 



» 



(n) 
2/5 


(n) 
— x 8 




= X 


(n) 


(n) 
— x 10 




= X 


(n) 
2/9 


(n) 
— S 12 




= X 



» 



(n) 

15 ' 
in) 
18 > 
(n) 

16 ' 
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the recurrence system ()A.3|) takes the following form 



( {n+i) 


= eef(A^f, 


(n- 

y\ 


HI) 






= (B{ n) f, 


(n- 


U) 


= e(A^r, 




= e x A^Bf\ 


(ri- 
Ve 




_ jWfiW 

— ^2 u 2 ' 


(n+1) 

y 7 


= eejiB^) 2 , 


(n- 


HI) 


= e 1 Bi n) ci n 


(n+l) 











where 



A 



B 



A. 



B 



C 



29 p y$ l) + 20 p y$ l) 
.(») 



(n) 
2/7 



V? + 2^ n) + 2y <"> + ##0 + 29 p y^ + ^ 



Noting that 



a2f„.( n )\2 _ „,("■)„, ( n ) 



(n) (n) 
(n) (n) 

y 7 % > 



we see that only five independent variables remain. 

It should be noted that if 0\ = 1, i.e. Ji = 0, then for the boundary condition a = r]i we have 



4"' 


— ^-2 




e 4 p B[ n \ 




- B in) 

— n 2 


— n 3 — 


c[ n) = 


Z 2 


_ y(jl) 

— Z 3 ' 







a 



(n) 



a 



in) 



so that 



z[ n+1) 

Z 3 



6(e$B[ n) ) 2 + 4(<%Bl n) )B[ n) + 2(9 + l)^) 2 



9(6*B{ n) ) 2 + Aie^B^B^ + 2(6 + 1)(B 



(n)No(n) 



,(n) )2 



Consequently, when #i = 1 for any boundary condition exists single limit Gibbs measure, namely, the 
unordered phase. So that the phase transition does not occur. 



Appendix B. Proof of the consistency condition 

In this section we show that the condition (|3.12fl and (j3.13j) are equivalent. Assume that (|3.12fl 
holds. Then inserting (|3.1Ujl into (|3.12|) we find 



nT II ^exp{/3Jicj(2;)(a(y)+a(z)) + /?J C j(y)a(z) 



(B.l) 



+h y a(y) + h z a(z)} = J| exp{h x a(x)}, 

xeWn-l 

» 



here given x G W n -\ we denoted = {y,z}, = {a(y),a(z)} and used = \J cr^K 

X&W n -l 
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Now fix x 6 W n -i and rewrite (jB.lj) for the cases a(x) = rji (i = 1,2) and a(x) = 773, and then 
taking their rations we find 

E CT M ={o . (2/)j0 . (z)} exp{/3Jir ?i (cr(y) + a{z)) + f3Ja(y)a{z) + h y a(y) + h z a(z)} 
E CT («) ={cr ( 2/ ) iCT(z )} exp{-pjir] 3 (a(y) + <r(*)) + 0Ja(y)a(z) + hycr{y) + h z a(z)} 
(B.2) = exp^}. 

Now by using l|2.9() from (|ij.2|) we get 
(B.3) A = F(h' y , h' z ), e h '*,* = F({h' y )\ [h! z f). 

From the equality HB,3|) we conclude that the function h = {h x = (h x i, h x ^) ■ x £ V} should satisfy 

(i3~m 

Note that the converse is also true, i.e. if (|3.13|) holds that measures defined by (|3,1U|) satisfy the 
consistency condition. Indeed, the equality (|3.13|) implies ()B.3|) . and hence (|B.2|) . From (|B.2|) we 
obtain 

exp{/?Jir/i(o-(y) + a(z)) + (3Ja(y)a(z) + h y a(y) + h z a(z)} 

ai n) ={a(y),a(z)} 

= a(x) exp{r]ih x }, 
where i = 1,2,3 and a(x) is some function. This equality implies 

Yl expipj^aiy) + a{z))a{x) + (3Jo{y)a{z) + h y a{y) + h z a(z)} 

(B.4) = Yl a{x)exp{a(x)h x }. 

Writing A n = Y\ xGW a{x) from (|B.4|) one gets 

(B.5) Z^^^^^K-^ZW^^K^V/'). 

<r(n) 

Taking into account that each > 1 is a probability measure, i.e. 

£ ^M (n) K-iv ff W) = 1, J2 M (n-1 Wi) = 1. 

from (|B.5j) we infer 

(B.6) Z^^An^ = Z< n », 

which means that (|3.12|) holds. 
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